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FINAL  REPORT  FOR  GRANT  AFOSR-87-0337 

High  resolution  SQUID  Magnetometry  for  Non-destructive  Evaluation 

John  P.  Wikswo,  Jr.,  Principal  Investigator 

Electromagnetics  Laboratory,  Department  of  Physics  and  Astronomy 
Vanderbilt  University,  Box  1807  Station  B,  Nashville,  TN  37235 

ABSTRACT 

Superconducting  QUantum  Interference  Device  (SQUID)  magnetometers  offer  promise  as  multi-mode 
instruments  capable  of  obtaining  high  resolution  images  of  extremely  low  frequency  injected  currents  or  eddy  currents, 
and  they  can  be  conrigured  to  image  the  magnetic  susceptibility  of  titanium,  aluminum  and  nonmetallic  composites. 
While  high  resolution  SQUID  magnetometers  will  generally  be  noisier  than  conventional  SQUIDs,  the  small  coils  and 
reduced  coil-to-source  spacing  more  than  compensate  to  provide  low-noise,  high-tesolution  images.  The  high  spatial 
resolution  which  can  be  obtained  with  SQUID  magnetometers,  the  unparalleled  sensitivity  of  SQUIDs  at  low 
frequencies,  the  ability  to  measure  weak  perturbations  in  strong  applied  magnetic  fields,  and  the  ability  to  discriminate 
against  external  sources  of  noise  should  make  SQUID  magnetometers  well  suited  for  NDE  of  deep  flaws  in  aluminum 
and  titanium  aerostructures.  To  explore  SQUID  NDE,  we  have  developed  research  facilities  that  include  the  high- 
resolution  MicroSQUID  magnetometer,  a  magnetic  shield,  a  scanning  stage,  and  a  computer-based  control  and  data 
acquisition  system.  Using  this  instrumentation,  we  have  imaged  magnetic  fields  produced  by  sources  as  varied  as 
intrinsic  currents  due  to  corrosion  or  Johnson  noise,  remanent  magnetization  from  ferromagnetic  contamination,  flaw- 
induced  perturbations  in  either  injected  current  or  eddy  currents  induced  by  an  AC  field,  Johnson  noise  currents  in 
a  copper  ring,  persistent  currents  in  high  transition-temperature  superconductors,  distributions  of  dia-  or  paramagnetism 
in  an  AC  or  DC  magnetic  field,  and  surface  flaws  decorated  with  a  paramagnetic  tracer.  In  support  of  the 
experimental  studies,  we  have  developed  analytical  and  numerical  models  for  the  simulation  of  flaws  with  several 
geometries  inside  thick  and  thin  current-carrying  plates  and  thin-walled  tubes,  and  have  demonstrated  that  two-dimen¬ 
sional  magnetic  images  can  be  deconvolved  into  images  of  current  or  magnetization  by  filtering  techniques,  finite 
element  models,  lead  field  analyses,  and  maximum  entropy  methods.  We  find  that  potential  applications  of  SQUID 
NDE  to  the  aging  aircraft  problem  include  high  resolution  imaging  of  flaws  and  corrosion  deep  inside  second  and  third 
layer  structures  and  turbine  rings,  and  magnetic  susceptibility  imaging  of  corrosion,  nonmetallic  composites  and 
magnetically-decorated  surface  flaws  in  ceramics.  However,  a  critical  assessment  of  the  capabilities  of  SQUID  NDE 
will  require  developing  quantitative  measures  of  SQUID  NDE,  such  as  the  probability  of  detection  of  various  sizes 
and  locations  of  flaws.  As  instruments  built  specifically  for  NDE  become  more  advanced,  and  as<«Kperience  is  gained 
worldwide  in  SQUID  NDE,  specific  applications  will  undoubtedly  be  found  for  which  the  capabilities  of  SQUIDs  are 
unmatched;  whether  or  not  this  is  sufficient  to  lead  to  widespread  utilization  of  SQUIDs  for  aircraft  NDE  will  then 
be  a  matter  of  economics,  ease  of  use,  and  perception  by  the  user  community. 

1.  INTRODUCTION 

While  SQUID  magnetometry  has  been  successfully  applied  to  a  broad  spectrum  of  problems  in  physics, 
biology,  engineering,  and  geology,  there  has  been  only  limited  use  of  SQUIDs  for  non-destructive  evaluation.  The 
first  demonstrations  were  by  Weinstock  and  Nisenoff'  in  1985.  Since  then,  SQUIDs  have  been  used  to  examine  or 
detect  ferromagnetic  plates^’^’^’*,  naval  mines®,  and  wrist  watches’;  as  an  amplifier  to  detect  eddy  current  signals®, 
for  monitoring  corrosion  in  buried  gas  pipes’  ’”;  and  for  fundamental  measurements  on  the  Barkhausen  effect"  ”. 
The  work  on  corrosion  on  pipes  is  now  l^ing  performed  using  fluxgate  magnetometers”.  There Jiave  been  only  two 
SQUID  NDE  review  articles”'”.  The  limited  effort  on  SQUID  NDE  has  occurred  for  several  reasons:  SQUID 
magnetometers  typically  have  more  sensitivity  than  is  necessary  for  most  NDE  applications;  SQUID  pickup  coils,  ordi¬ 
narily  optimized  for  clinical  biomedical  measurements  on  humans,  are  too  large  and  too  far  from  the  sample  to  provide 
the  requisite  spatial  resolution  to  compete  with  state-of-the-art  eddy  current  and  ultrasound  techniques;  and  few  of  the 
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research  groups  doing  SQUID  NDE  have  programs  that  simultaneously  develop  instrumentation  and  mathematical 
models  while  conducting  experiments  that  push  the  understanding  of  both  the  magnetic  fields  and  their  measurement. 

One  of  the  unfortunate  misperceptions  of  the  NDE  user  community  is  that  SQUID  magnetometers  are  suitable 
only  for  laboratory  analysis.  In  part,  this  error  arises  from  the  requirement  of  SQUID  magnetometers  for  liquid  heli¬ 
um.  There  are  several  notable  examples  that  demonstrate  the  success  of  SQUIDs  outside  the  laboratory.  For  many 
years,  SQUID-based  rock  magnetometers  have  even  been  installed  in  deep  sea  drilling  ships,  and,  through  the  use  of 
a  nitrogen-temperature  cryocooler,  have  a  liquid  helium  hold-time  of  almost  a  year"'’.  A  compact  SQUID  magne¬ 
tometer  system,  recently  available  commercially,  is  capable  of  being  operated  at  the  end  of  a  robot  manipulator  arm 
in  a  magnetically-harsh  industrial  environment’’.  A  70  gigahertz  sampling  oscilloscope  with  less  than  a  one-minute 
cool  down  time  has  been  sold  commercially,  and  a  prototype  version  has  been  produced  using  a  closed-cycle  refrig¬ 
erator'*.  In  fairness  to  the  cryogenically  inexperienced,  the  need  to  transfer  liquid  helium,  its  cost,  and  the  logistical 
problems  with  securing  delivery  of  liquid  helium  to  remote  sites  have  all  contributed  to  the  slow  acceptance  of 
SQUID-based  devices  outside  the  laboratory. 

The  discovery  of  high  temperature  superconductors  with  transition  temperature  above  100  K  will  change  this 
situation  completely.  Steady  progress  has  been  made  towards  the  fabrication  of  high  temperature  superconducting 
components  and  devices  by  a  number  of  companies,  including  DuPont,  Conductus,  TRW,  IBM,  Biomagnetic 
Technologies,  and  American  Superconductor,  as  well  as  several  university  laboratories  and  government  agencies  such 
as  the  National  Institute  of  Standards  and  Technology.  Several  laboratories  have  produced  DC  SQUID  devices  with 
lower  noise  than  helium-temperature  RF  SQUID  magnetometers,  and  there  are  predictions  that  the  high-frequency 
sensitivity  of  high-temperature  SQUIDs  will  soon  approach  that  of  commercially-available  low  temperature  SQUIDs. 
A  high-Tj  SQUID,  Mr.SQUID“",  is  now  available,  as  an  off-the-shelf  product,  for  educational  use”,  and  by  year’s-end 
a  version  of  this  device  with  complete  feedback  electronics  may  be  available  as  a  simple  yet  complete  SQUID  mag¬ 
netometer  system  whose  sensitivity  will  far  exceed  that  of  commercial  fluxgate  magnetometers.  Thus  in  the  course 
of  the  next  several  years,  it  is  reasonable  to  assume  that  a  spectrum  of  high-T^  SQUID  magnetometer  systems  will  be 
available  commercially,  if  not  on  a  production  basis  definitely  on  a  custom-fabrication  one. 

The  discovery  of  high-temperature  superconductors  has  also  initiated  major,  new  research  programs  in 
cryogenic  technology,  since  100  K  can  be  readily  achieved  with  an  inexpensive,  portable  refrigerator.  It  would  not 
be  surprising  to  find  a  SQUID  magnetometer  mounted  on  a  miniature  refrigerator  within  3  to  4  years,  or  even  sooner. 
While  it  is  unlikely  that  high-T^  SQUIDs  and  their  associated  refrigerators  will  provide  the  absolute  sensitivity  pres¬ 
ently  afforded  by  the  low-T^  SQUID  systems  being  manufactured  for  magnetoencephalography,  it  should  also  be  clear 
from  the  data  presented  in  this  paper  that  most  SQUID  NDE  signals  are  vastly  larger  than  biomagnetic  ones.  As  an 
example  of  readily  achieved  sensitivity,  a  high-T<,  SQUID  with  a  2  mm’  area  hole  in  the  SQUID  washer  will  have  a 
flux  noise  of  10“*<|)yV’Hz,  which  would  translate  to  a  field  noise  of  10  '^  T// Hz,  equivalent  to  that  presently  achieved 
by  the  liquid  helium-cooled  MicroSQUID  magnetometer  in  use  at  Vanderbilt’”. 

The  key  concern,  of  course,  is  not  only  SQUID  noise,  but  environmental  noise.  In  the  past  several  years,  there 
have  been  some  noteworthy  accomplishments  in  gradiometer  design  and  active  noise  cancellation.  Because  the  flaws 
being  examined  in  flight-line  NDE  of  aircraft  would  be  at  most  several  centimeters  beneath  the  outer  surface, 
gradiometer  baselines  can  be  of  comparable  dimension,  which  is  substantially  shorter  than  that  used  in  many  bio¬ 
medical  systems.  The  addition  of  reference  magnetometer  channels  to  sense  ambient  noise  and  the  use  of  computer- 
controlled  AC  cancellation  of  drive  current  crosstalk  has  made  it  possible  to  operate  SQUID  NDE  systems  in  a  particu- 
larly-hostile  industrial  environment”.  In  the  case  of  current  injection  or  eddy  currents,  it  is  simply  possible  to  increase 
the  drive  current  into  the  sample  until  the  signal  rises  from  the  noise,  within  the  limitations  imposed  by  SQUID 
dynamic  range;  all  of  the  Vanderbilt  data  presented  later  in  this  paper  utilized  total  drive  currents  of  less  than  100 
milliamps!  The  Magneto-Optic  eddy  current  Imager  (MOI)”  demonstrates  that  magnetic  imaging  is  feasible  on  the 
flight  line;  SQUID  magnetometry  may  offer  the  potential  to  achieve  flaw  detection  with  even  greater  reliability  and 
sensitivity. 
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All  of  these  factors,  i.e.  high  temperature  superconductivity,  cryogenic  refrigeration,  the  relative  modest 
sensitivity  requirements  of  SQUID  NDE,  and  the  demonstrated  capabilities  of  noise  cancellation  techniques  combine 
to  make  flight-line  SQUID  NDE  feasible.  However,  preliminary  aircraft  NDE  studies  to  date  only  show  that  the 
technique  has  promise;  definitive  demonstration  of  the  sensitivity  and  reliability  of  SQUID  for  aircraft  has  yet  to  be 
accomplished.  As  we  will  describe  below,  the  work  at  Vanderbilt  has  been  directed  towards  the  continued 
development  of  low  temperature  instrumentation  and  techniques  to  provide  a  rigorous  evaluation  of  the  potential 
strengths  and  weaknesses  of  SQUID  NDE.  The  subsequent  development  of  specialized  SQUID  NDE  instruments 
suitable  for  use  outside  the  laboratory,  possibly  utilizing  high-T^  instruments,  should  allow  the  transfer  of  this  new 
technology  to  the  flight-line. 


2.  THE  VANDERBILT  SQUID  NDE  EFFORT 

The  Vanderbilt  group  has  examined  both  experimentally  and  theoretically  a  variety  of  problems  in  SQUID 
NDE.  We  have  developed  a  specialized  Magnetic  Imaging  facility  for  studying  the  magnetic  fields  from  biological 
and  non-biological  systems.  We  have  imaged  the  magnetic  fields  produced  by  a  wide  variety  of  sources.  In  this 
review,  we  will  discuss  our  various  NDE  measurements  in  the  context  of  aircraft.  The  parallel  effort  in  biophysics, 
particularly  the  imaging  of  cardiac  action  currents,  is  discussed  elsewhere^'^.  We  have  also  developed  analytical 
and  numerical  models  that  describe  our  measured  magnetic  fields,  and  we  have  developed  inverse'algorithms  to  allow 
us  to  deconvolve  the  fields  to  obtain  images  of  the  sources,  whether  they  be  distributions  of  current  density  or 
magnetization.  We  now  describe  each  of  these  aspects  of  our  work. 

2.1  Instrumentation 

The  key  hardware  systems  in  the  Vanderbilt  Magnetic  Imaging  Facility  are  the  MicroSQUID  magnetometer^, 
a  magnetic  shield^,  a  scanning  staged,  and  a  computer  control  and  data  acquisition  system^.  The  MicroSQUID 
magnetometer  in  Fig.  1  overcomes  the  limitations  of  previous  biomedical  SQUID  systems  by  utilizing  small  (3  mm) 
diameter  coils  located  close  (within  1 .5  mm)  to  room  temperature  samples.  \^ile  the  90-l(X)  ff/Hz'^  noise  in  Micro¬ 
SQUID  is  larger  than  in  conventional  SQUID  magnetometers,  the  increased  resolution  and  the  reduced  coil-to-source 


Figure  1.  The  VanderbUt  MicroSQUID  magnetometer.  Two  views  of  the  tail  of  the  Dewar,  showing  the  close 
spacing  of  the  pickup  coils  to  the  outside  of  the  dewar“. 
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Figure  2.  A  test  of  the  FTT  imag* 
ing  algorithm,  a)  A  printed  circuit 
pattern  carrying  a  100  fiA  DC  cur¬ 
rent.  b)  The  magnetic  field  map 
measured  by  MicroSQUID  with  the 
3  mm  diameter  pickup  coils  posi¬ 
tioned  2.5  mm  above  the  circuit,  c) 
The  current  image  obtained  from  the 
data  in  b)  using  unconstrained  FFT 
deconvolution  (12.5  /tA  per  line). 
Ads^ted  from  Ref.  25. 


spacing  more  than  compensate  to  provide  low-noise,  high-resolution  images^.  Multiple  SQUHT&iannels  allow  real¬ 
time  spatial  filtering  that  is  optimized  for  the  particular  source  being  studied,  while  at  the  same  time  iti^roving  the 
signal-to-noise  ratio.  The  magnetic  shield  provides  sufficiently  low  noise  (8  fT  at  10  Hz)  and  high  60  Hz  attenuation 
(10^)  that  it  serves  as  an  essentially  noise-free  SQUID  NDE  test  facility.  The  first  computerized,  nonmagnetic  transla¬ 
tion  stage  for  scanning  samples  beneath  MicroSQUID,  with  its  50  x  50  mm  range  and  1  mm/s  velocity,  has  just  been 
replaced  by  a  high  spe^  (20  cm/s),  22  cm  x  66  cm  stage. 

2.2  Image  processing 

The  success  and  accqjtance  of  many  NDE  techniques  such  as  radiography,  ultrasound  C-scans,  ultrasonic 
microscopy,  and  X-ray  tomography  have  been  aided  by  the  presentation  of  two-or  even  three-dimensional  images  of 
the  test  object  and  the  flaws  that  it  contains.  However,  Ae  lack  of  an  image  has  not  prevented  the  widespread 
acceptance  of  conventional  eddy-current  and  ultrasound  instruments,  for  which  the  nature  of  the  measurement  often 
provides  a  characteristic,  qualitative  signal  or  signature  rather  than  a  quantitative  image  of  thcTlaw.  With  SQUID 
NDE,  the  uniformity  of  the  iqiplied  magnetic  field  or  current,  the  accuracy  with  which  the  SQUID  can  nuq)  the  field 
above  the  test  object,  and  the  ease  with  which  any  distortion  of  these  fields  by  the  unflawed  portions  of  the  object 
can  be  eliminated  make  quantitative  flaw  imaging  an  attractive  capability  of  the  technique. 
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We  have  previously  demonstrated  that  two-dimensional  magnetic  images  can  be  deconvolved  into  images  of 
current  or  magnetization  by  filtering  techniques  carried  out  in  the  spatial  frequency  domain^,  as  was  done  for  the 
test  pattern  in  Fig.  2.  However,  with  this  approach,  it  is  difficult  to  utilize  a  priori  information  such  as  the  known 
geometry  of  the  irregular  borders  of  the  sample.  The  need  to  constrain  images  can  be  readily  met  through  a  finite 
element  imaging  technique^-”,  illustrated  in  Fig.  3.  Other  imaging  approaches  include  lead  field  analyses^, 
maximum  entropy  methods,  iterative  approaches”,  and  others  that  are  under  investigation  at  Vanderbilt  and 
elsewhere. 


Different  deconvolution  algorithms  have  their  own  merits  in  dealing  with  certain  kinds  of  sources.  As  we 
will  show  in  the  following  sections,  we  have  applied  appropriate  algorithms  to  magnetic  field  data  produced  by  a 
planar  current-carrying  conductor  with  flaws,  and  by  the  diamagnetic  magnetization  of  plexiglass.  As  an  example 
of  the  differences  between  techniques,  the  finite  element  method  can  constrain  the  solution  to  a  finite  extent,  but  is 


Figure  3.  FFT  and  finite  element  deconvolution,  a)  A  numerical  test  current  pattern,  b)  An  unconstrained  FFT 
reconstruction  of  the  currents  from  simulated  magnetic  fields  3  mm  above  the  pattern  with  20%  noise  (Mean  square 
deviation  =  0.42).  c)  A  constrained  finite  element  mesh  defining  the  conductor  borders,  d)  A  finite  element  recon¬ 
struction  from  the  same  magnetic  data  (MSD  =  0.122).  Adapted  from  Ref.  29. 
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sensitive  to  the  imperfections  in  the  input  data.  The  maximum  entropy  method  is  more  stable  for  fitting  noisy  data 
but  may  not  find  the  best  solution  for  data  with  a  high  signal-to-noise  ratio.  For  cylindrical  geometries,  a  lead  field 
analysis  provides  a  simpler  algorithm.  Each  of  these  algorithms  has  advantages  for  particular  problems  while  being 
inadequate  for  others.  However,  with  all  the  approaches,  the  deconvolved  source  images  are  easier  to  interpret  than 
the  original  magnetic  field  data^’.  Whether  a  complete  image  or  simply  a  characteristic  signal  are  required  to  detect 
a  flaw  may  depend  upon  the  specific  application  being  addressed;  in  either  case,  our  preliminary  studies  indicate  that 
one  of  the  advantages  of  SQUID  NDE  is  that  it  provides  quantitative  data  that  can  be  viewed  directly  or  analyzed  by 
measurement  models. 


5  nT  contours 


1  nT  contours 


Figure  4.  Cancellation  of  background  signals,  a)  A  25  mm  x  150  mm  x  32  /xm  copper  sheet  with  a  3  mm 
diameter  hole  at  its  center.  A  7.5  mA,  1.6  Hz  current  was  applied  through  contacts  at  each  end.  The  dotted  square 
represents  the  region  scanned  by  MicroSQUID.  b)  An  isofield  plot  of  the  normal  component  of  the  magnetic  field 
over  the  region  indicated  by  the  dotted  square  in  (a).  The  major  variation  was  due  to  the  field  from  the  plate  edges, 
while  the  small  variation  at  the  center  was  due  to  the  hole,  c)  The  canceling  technique  as  applied  to  the  conductor 
in  (a).  The  unflawed  canceling  plate  was  450  /xm  beneath  the  conductor  in  (a),  and  carried  an  equal  and  opposite 
current.  The  two  conductors  were  coimected  to  the  power  supply  by  a  coaxial  cable,  d)  An  isofield  plot  for  the 
magnetic  field  in  the  dotted  square  in  (c).  Note  the  increased  scale  as  compared  to  (b).  Adapted  from  Ref  32. 
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2.3  Experimental  studies 

2.3.1  Injected  current  imaging.  Weinstock'  first  demonstrated  the  ability  of  a  SQUID  magnetometer  to  detect 
flaws  in  conducting  objects  by  observing  the  magnetic  fields  associated  with  the  perturbation  in  an  applied  electric 
current  as  it  is  deflected  by  the  flaw.  Unfortunately,  when  this  technique  is  applied  to  flaws  in  planar  conductors, 
the  signal  from  the  flaw  can  be  masked  by  that  of  the  currents  in  the  conductor  and  the  lead-in  wires,  even  though 
the  flaw  signature  is  well  above  the  noise^'.  This  limitation  has  been  overcome  by  a  technique  developed  at  Vand¬ 
erbilt  in  which  an  unflawed  conductor  carrying  an  equal  and  opposite  current  distribution  is  placed  adjacent  to  the 
test  surface^^.  As  shown  in  Fig.  4,  this  allows  the  current  to  be  delivered  by  a  magnetically-silent  coaxial  cable, 
and  eliminates  most  of  the  background  magnetic  field  due  to  the  currents  in  the  unflawed  sections  of  the  test  plate. 


5  mm 


5  mm 


X,  mm 

Figure  5.  SQUID  imaging  of  hidden  corrosion,  a)  A  simulated  lap  joint  in  an  aircraft  wing  surface,  consisting 
of  two  insulated  aluminum  plates.  Standard  eddy  current  analysis  indicates  a  2%  mass  loss  in  the  corroded  area. 
A  10  Hz  58  mA  AC  current  was  injected  into  the  plate,  b)  and  c)  Magnetic  images;  the  dashed  line  shows  the  hidden 
flaw.  From  Ref.  34. 
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Figure  6.  SQUID  detection  of  second-layer  cracks,  a)  A  test  sample  of  two  7075-T6  alCufRiium  panels  bolted 
together  by  six  aluminum  pins  and  nuts.  Each  panel  is  250  x  100  x  2.4  mm  thick.  For  three  pins,  second-layer 
cracks  are  simulated  by  EDM  slots  in  the  bottom  layer  adjacent  to  the  rivets.  The  slot  width  was  less  than  0.4  mm. 
The  SQUID  pick-up  coil  is  located  3  mm  above  the  rivets  (4.5  mm  above  the  plate).  47  Hz  30  mA  current  was 
injected  into  die  plate  and  the  magnetic  field  3  mm  above  each  pin  was  measured  with  MicroSQUID.  c),  d),  f),  and 
g)  Isofield  contour  plots  for  four  different  rivets,  b)  and  e)  Corresponding  field  profiles  through  the  rivets.  The 
crack  in  (c)  is  12  mm  (1/2  in)  long  and  transverse  to  the  cunent;  while  in  (f)  it  is  parallel  to  the  current.  Scans  (d) 
and  (g)  are  for  adjacent,  unflawed  holes  with  rivets.  The  differences  between  the  flawed  and  unflawed  holes  with 
rivets  are  evident  in  the  current  profiles  (b)  and  (e)  along  a  horizontal  line  through  the  holes.  Adapted  from  Ref.  34. 


In  collaboration  with  Lxickheed’^,  we  have  begun  to  apply  our  techniques  to  the  aging  aircraft  problem. 
Figure  5  shows  magnetic  data  recorded  from  Lxickheed-prepared  aluminum  samples  that  simulate  lap  joints  in  aircraft 
wings^.  Injection  of  a  58  mA,  10  Hz  current  allowed  detection  of  a  hidden,  16  mm  radius  corrosion  area,  whose 
mass  loss  was  approximately  2%  according  to  eddy  current  analysis.  As  shown  in  Fig.  6,  our  technique  is  also  capa¬ 
ble  of  discriminating  between  the  current  perturbation  caused  by  an  aluminum  fastener  and  that  caused  by  a  6.4  nun 
long  EDM  slot  immediately  adjacent  to  the  fastener.  By  varying  the  frequency  of  the  injected  current,  it  is  possible 
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to  adjust  the  maximum  depth  sensitivity  of  the  technique:  at  10  Hz,  the  skin  depth,  5,  is  4  cm  in  7075  Al,  and  so 
sensitivity  is  limited  only  by  the  1/r^  dependence  of  the  field  due  to  the  flaw.  As  we  will  discuss  later,  the  width  of 
the  signature  is  directly  related  to  the  depth  of  the  flaw,  and  the  strength  of  the  signature  to  its  size  and  orientation. 
The  deconvolution  algorithms  discussed  above  can  also  be  used  to  obtain  images  of  the  current  distribution  in  the 
sheet  for  both  flawed  and  unflawed  second-layer  rivets,  presented  in  Fig.  7,  which  show  how  the  current  pattern  is 
both  displaced  to  the  left  and  covers  a  larger  region  because  of  the  EDM  slot  on  the  left  side  of  the  rivet. 
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Figure  7.  Field  maps  and  current  images  of  rivets.  Measured  magnetic  field  contour  map  (above)  and  computed 
current  density  arrow  map  (below)  due  to  current  injection  into  a  pair  of  conducting  plates  containing  rivets  without 
flaws  (left)  and  with  second-layer  flaws  (right).  From  Ref.  28. 
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We  have  inined  a  1.1  mm  thick  test  plate^*  with  a  1  mm  long  crack  located  entirely  within  the 
countersunk  region  for  a  flush-mounted  rivet^*.  An  injected  current  density  of  only  2.5  mA/cm^  was  required  to 
produce  the  images  in  Fig.  8.  By  imaging  the  magnetic  field  for  current  parallel  to  the  crack,  we  can  estimate  the 
signature  from  an  unflawed  rivet  hole.  There  is  a  clear  difference  between  the  magnetic  images  8d  and  8e,  with  the 
peak  to  peak  amplitude  and  the  separation  of  the  peaks  in  8f  being  significantly  larger  for  the  case  where  the  current 
has  to  flow  around  the  crack.  The  pair  of  loops  in  the  FFT  current  image  in  Fig.  8g  may  be  the  actual  image  of  the 
two  cracks,  but  the  more  detailed  analysis  of  the  imaging  process  will  be  required  before  we  are  able  to  interpret  this 
image  fully.  Analysis  of  images  recorded  at  a  frequency  of  several  hundred  Hz  suggests  that  it  may  be  possible  to 
achieve  selective  depth  sensitivity  that  enhances  the  crack  signature^’,  but,  again,  more  theoretical  analysis  is  required. 


X,  mm  X,  mm  x,  mm 


J  =  2.54  mA/cm^  DC  Current  Injected 

Figure  8.  SQUID  detection  of  cracks  within  a  rivet  countersink.  Test  specimen  with  EDM  machined  slots  that 
do  not  extend  beyond  the  countersink  of  the  rivet  hole^’.  a)  and  b)  Two  orientations  of  applied  current,  c)  Sample 
cross-section,  d)  and  e)  Isofield  maps  for  currents  injected  perpendicular  and  parallel  to  the  slot,  f)  Field  profiles 
showing  the  detectability  of  the  slots,  g)  and  h)  Current  images  computed  from  (d)  and  (e)  using  the  FFT  algorithm, 
i)  Current  profiles  through  the  flaw. 
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Figure  9.  SQUID  imagiiig  of  deep  flaws  in  metals,  a)  A  4.7  mm  radius  spherical  cavity  inside  a  305  x  88  x  25.4 
mm  rectangular  brass  bar  carrying  70  mA  of  DC  current  along  its  length.  The  cavity  is  located  7.9  mm  below  the 
top  surface,  b)  Surface  and  c)  isofield  contour  maps  of  the  magnetic  field.  d,e,f)  Maps  for  10,  130,  and  200  Hz  AC 
currents,  respectively,  g)  Theoretical  and  experimental  field  profiles.  Adapted  from  Ref  37. 
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Subsurface  flaws  in  electrically  conducting  components  are  difficult  to  detect  by  conventional  eddy  current 
techniques  because  of  the  skin  depth,  6,  at  high  frequencies.  However,  SQUID  magnetometers  are  capable  of  imag¬ 
ing  DC  and  low-frequency  magnetic  fields,  thereby  allowing  sources  to  be  detected  deep  inside  conducting  samples. 
A  spherical  flaw  beneath  the  surface  of  a  brass  bar  has  been  detected  by  injecting  both  DC  and  low-frequency  AC 
current.  Figure  9  shows  magnetic  images  for  a  spherical  cavity  with  a  4.7  mm  radius  inside  a  305  mm  long  x  88 
mm  wide  X  25.4  mm  thick  rectangular  brass  bar^’.  At  DC,  the  flaw  is  readily  imaged.  By  injecting  AC  current 
with  several  frequencies,  we  can  estimate  the  depth  of  the  flaw.  At  frequencies  above  200  Hz,  the  flaw  is  practically 
invisible,  indicating  the  potential  benefits  of  SQUID-based  extra  low  frequency  eddy  current  (ELF-EC)  techniques. 

2.3.2  Eddy  current  imaging.  While  the  injected-current  technique  is  useful  for  high-precision  measurements 
on  test  samples,  ihe  need  to  make  good  electrical  contact  with  the  sample  would  make  it  difficult  to  use  the  technique 
on  actual  aircraft  stmctures,  particularly  if  they  are  painted.  We  have  adapted  a  standard  eddy-current  technique, 
also  used  in  the  Magneto-Optic  Imager^*-^’,  in  which  the  AC  magnetic  field  is  applied  by  a  sheet  conductor  parallel 
to  the  test  surface.  This  induces  a  large-extent  sheet  current  in  the  test  specimen,  and  thus  produces  flaw  perturbation 
fields  quite  similar  to  those  obtained  with  direct  current  injection^.  SQUID  data  from  this  technique  are  shown  in 
Figs.  10,  11  and  12,  confirming  that  this  approach  is  valid  for  hidden  corrosion,  closed  fatigue  cracks,  and  second- 
layer  cracks  at  rivets. 


Figure  10.  SQUID  eddy  current  image  of  corrosion,  a)  An  insulated,  2.5  mm  thick  aluminum  plate.  The  26  mm 
X  21  mm  corrosion  area  is  open  on  the  surface,  and  is  less  than  1  mm  deep.  Another  flat,  current-carrying  plate 
between  the  SQUID  and  the  test  sample  is  used  to  induce  the  eddy  currents,  whose  magnetic  field  is  shown  in  b). 
The  flaw  could  be  imaged  from  either  side. 
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Figure  11.  SQUID  imaging  of  a  dosed  fatigue  crack.  A  MicroSQUID  image  of  the  magnetic  field  produced  by 
the  perturbation  in  induced  currents  that  results  from  a  closed  fatigue  crack  half-way  across  a  25  mm  x  2.3  mm 
aluminum  strip  with  a  rivet  hole”,  a)  Plan  and  b)  cross-section  views,  c)  Isofield  contour  map. 

2.3.3  Magnetic  susceptibility  imagine.  The  SQUID  NDE  techniques  that  we  have  described  so  far  are  limited 
to  conducting  objects.  While  these  techniques  should  be  well-suited  for  the  detection  and  quantification  of  second- 
layer  cracks  and  corrosion  in  aging  aluminum-skin  aircraft,  they  will  be  of  little  use  for  aircraft  fabricated  using 
advanced  non-metallic  composites  or  ceramic  materials.  Ultrasound  C-scans  are  useful  for  finding  internal  voids  in 
planar  sheets  of  fiber/resin  composites,  but  the  utility  of  ultrasound  as  a  NDE  tool  is  compromised  for  complex 
structured  shapes,  largely  because  of  the  difficulty  in  maintaining  the  transducer  in  good  acoustical  contact  with  the 
material,  because  of  internal  reflections  and  scattering  of  the  acoustic  waves  by  curved  surfaces  and  internal 
boundaries,  and  because  of  lack  of  contrast  between  some  inqiurities  and  the  base  material.  X-ray  techniques  are 
limited  by  the  low  x-ray  attenuation  coefficient  of  many  of  these  materials,  resulting  in  minimal  contrast  between  a 
void  and  the  adjacent  solid  material.  Other  studies  demonstrated  that  SQUID  mapping  of  remanent  magnetic  fields 
can  be  used  to  study  objects  with  intrinsic  ferromagnetic  impurities  or  extrinsic  ferromagnetic  decoration^'-^-*', 
but  the  inverse  imaging  problem  for  these  measurements  is  insolvable  because  of  the  lack  of  knowledge  regarding 
the  direction  of  magnetization  of  the  ferromagnetic  material. 

Several  of  these  difficulties  may  be  overcome  using  magnetic  susceptibility  imaging  and  tomography 
techniques  developed  at  Vanderbilt^^'^^’*’”’“.  We  have  shown  experimentally  and  theoretically  how  SQUID 
magnetometers  can  be  used  in  conjunction  with  an  applied  magnetic  field  to  image  spatial  variations  in  the  magnetic 
susceptibility  of  diamagnetic  and  paramagnetic  objects.  In  order  to  demonstrate  the  feasibility  of  this  technique,  we 
have  attached  a  pair  of  Helmholtz  coils  to  our  high  resolution  MicroSQUID  magnetometer  so  that  we  can  apply  a 
nearly  uniform  field  of  0.3  mT,  2  mm  below  the  3  mm  diameter  sensing  coils^.  By  placing  a  small  sample  in  this 
applied  field  and  scanning,  MicroSQUID  measures  the  local  magnetic  field  variations  caused  by  the  spatial  distribution 
of  magnetic  susceptibility  in  the  sample.  The  system  has  sufficient  sensitivity  to  image  susceptibility-induced  field 
changes  equal  to  3  x  lO  *'’  of  the  applied  field,  in  a  1  Hz  bandwidth. 
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Figure  12.  SQUID  eddy-current  detection  of  second  layer  cracks.  For  the  same  sample  as  in  Fig.  6,  a  flat,  cur¬ 
rent-carrying  plate  between  the  SQUID  and  the  test  sample  was  used  to  induce  340  Hz  eddy  currents  across  the 
sample  strip,  c),  d),  f),  and  g)  Isofield  contours  for  four  different  rivets,  b)  and  e)  Conesponding  field  profiles 
through  the  rivets.  The  crack  in  (f)  is  transverse  to  the  current  while  that  in  (c)  is  parallel  to  it.  The  holes  in  (d) 
and  (g)  are  unflawed.  The  differences  between  the  flawed  and  unflawed  holes  are  evident  in  the  current  profiles  (b) 
and  (e)  along  a  vertical  line  through  the  holes.  Adapted  from  Ref.  34. 
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To  test  the  capabilities  of  this  technique  for  imaging  two-dimensional  distributions  of  magnetic  susceptibility 
in  thin  samples,  we  imaged  complex  shapes  of  plexiglass,  scanned  in  a  uniform  applied  field.  Deconvolution  of  these 
data“  yielded  susceptibility  distributions  that  were  accurate  images  of  the  samples,  as  shown  in  Fig.  13.  We  found 
that  the  susceptibility  of  this  plexiglass  sample  was  -9.0  x  10'*  (SI)  and  demonstrated"”  that  our  system  is  sensitive 
to  susceptibility  contrasts  as  small  as  5  x  10  ’  (SI),  with  a  spatial  resolution  on  the  order  of  1  mm.  We  have  also 
shown  that  this  technique  can  be  used  to  obtain  high  quality  images  of  the  heterogeneous  magnetic  susceptibility  and 
remanent  magnetization  of  a  50  fim  thick  slice  of  volcanic  rock^.  The  performance  of  the  MicroSQUID  implemen¬ 
tation  of  this  concept  is  limited  by  the  strength  of  the  applied  field,  by  the  size  and  distance  of  the  sensing  coils  from 
the  sample,  and  by  SQUID  noise.  We  estimate  that  we  could  increase  spatial  resolution  by  a  factor  of  two,  and 
sensitivity  to  susceptibility  contrast  by  at  least  one  order  of  magnitude  by  using  an  imaging  susceptometer  whose 
SQUIDs  are  integrally  mounted  within  a  superconducting  magnet.  Such  a  system  in  presently  under  construction  for 
our  studies. 
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Figure  13.  A  magnetic  image  of  plexiglas.  The  sample  is  a  25.4  x  25.4  x  6.3  mm  high  block  of  plexiglas  (sus¬ 
ceptibility  =  -9.0  X  10"*,  SI)  with  five  1.8  mm  diameter  holes  drilled  through  it.  a)  The  magnetic  field  recorded 
above  the  specimen  with  no  applied  field,  b)  The  magnetic  field  measured  with  an  applied  fieW-of  114  /xT.  c)  The 
2-D  magnetization  image  calculated  by  FFT  deconvolution  of  the  magnetic  field  perturbations  measured  2.5  mm 
above  the  block  in  the  presence  of  an  applied  field.  The  five  peaks  do  not  extend  to  the  zero  magnetization  level 
because  of  the  FFT  algorithm  and  the  limiting  resolution  of  the  3  mm  sensing  coils  of  the  SQUID.  From  Ref.  44. 


We  have  also  shown  theoretically  that  the  application  of  the  magnetic  field  from  multiple  directions  and  the 
recording  of  the  diamagnetic  perturbation  signals  over  the  entire  surface  of  the  object  can  be  utilized  to  obtain  a  three- 
dimensional  tomographic  reconstruction  of  the  internal  susceptibility  distribution  of  the  object^’,  overcoming  the 
famous  non-uniqueness  of  the  magnetic  inverse  problem.  To  obtain  the  three-dimensional  susceptibility  distribution 
X,  we  solve  the  system  of  equations  B  =  G  *  X,  in  which  B  is  the  set  of  magnetic  field  measurements  and  G  is  the 
matrix  of  Green’s  functions  that  must  be  inverted.  In  many  cases,  the  solution  is  non-unique,  but  we  overcome  this 
problem  by  measuring  the  magnetic  field  in  several  different  planes  close  to  the  sample  while  applying  the  magnetic 
field  sequentially  from  multiple  directions.  The  resulting  system  of  equations  is  overdetermined  and  singular  value 
decomposition  techniques,  among  others,  can  be  used  to  reconstruct  the  3-D  susceptibility  distribution.  Our  prelimi¬ 
nary  studies  show  that  with  two  of  the  three  strategies  investigated,  the  algorithm  can  correctly  discriminate  between 
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voxels  with  susceptibilities  differing  by  only  5  x  10  *  (SI).  Since  this  preliminary  numerical  simulation  utilized  only 
a  small  number  of  sample  orientations  and  a  single  magnet/SQUlD  configuration,  in  contrast  to  the  hundred  or  more 
projections  for  X-ray  tomography,  we  anticipate  a  significant  improvement  in  both  the  sensitivity  and  spatial 
resolution  of  magnetic  susceptibility  tomography  when  more  advanced  algorithms  and  measurement  protocols  are 
developed. 

Detection  of  surface  flaws  in  critical  metallic  components  is  currently  performed  by  eddy  current  techniques, 
or  by  visual  analysis  of  fluoro-penetrants.  Acoustic  microscopy  is  used  as  a  research  tool.  However,  these  methods 
are  often  limited  by  low  sensitivity,  poor  spatial  resolution,  or  slow  throughput.  Some  are  ill-suited  for  ceramics. 
All  of  these  limitations  are  addressed  with  an  alternative  technique  involving  decoration  of  the  non-ferromagnetic  sam¬ 
ples  with  a  ferromagnetic  or  superparamagnetic  penetrant  and  subsequent  scanning  with  a  high-resolution  mag¬ 
netometer  or  susceptometer.  Figure  14  presents  our  preliminary  data  that  demonstrated  that  this  technique  is  sensitive 
enough  to  detect,  with  a  10  to  1  signal-to-noise  ratio,  defects  at  least  as  small  as  250  /xm  long  x  50  fim  wide  x  75 
fim  deep  in  a  nickel  test  block"*^.  Applied  field  strength  and  SQUID  noise  will  ultimately  limit  sensitivity,  and  sensing 
coil  characteristics  will  limit  spatial  resolution.  We  expect  that  the  imaging  susceptometer  presently  under  construc¬ 
tion  will  improve  the  sensitivity  to  the  point  where  a  large-coil  SQUID  may  provide  both  rapid  and  sensitive  scanning 
of  ceramic  bearings  in  which  microscopic  surface  flaws  have  been  decorated  with  superparamagnetic  tracers. 


Figure  14.  Surface  decoration  with  superparamagnetic  microspheres.  Two  plots  of  the  magnetic  field  measured 
2.5  mm  above  a  block  of  nickel  alloy  containing  three  EDM  surface  defects  that  are  5(X)  ^ni  long  and  50  fim  wide 
with  various  depths  (left;  150  fim;  middle:  230  fim;  right  300  fim),  that  is  exposed  to  a  171  applied  field.  A  solu¬ 
tion  of  superparamagnetic  microspheres  was  applied  uniformly  to  the  block  and  its  surface  was  then  wiped  clean. 
The  field  from  the  susceptibility  of  the  block  was  removed  by  subtracting  the  field  measured  before  application  of 
the  spheres  from  that  measured  afterwards.  From  the  strength  of  the  signal,  the  small  defect  in  the  upper  left  corner 
was  found  to  have  a  volume  of  1.9  ±  0.6  x  10  '^  m^  and  was  subsequently  identified  through  microscopic 
examination  to  be  a  relatively  deep  scratch  approximately  as  long  as  the  left-most  EDM  slot.  Adapted  from  Ref.  45. 


Magnetic  susceptibility  imaging  offers  promise  as  a  new  method  for  imaging  composites,  either  through  sur¬ 
face  or  three-dimensional  tomographic  imaging  to  determine  regional  variations  in  susceptibility  arising  from  voids, 
unpolymerized  resin,  or  differences  in  the  relative  concentrations  of  chemical  components,  or  through  the  use  of 
ferromagnetic  or  superparamagnetic  penetrants.  Since  the  approach  is  relatively  new  and  untested  in  realistic  aircraft 
specimens,  it  is  premature  to  speculate  whether  this  approach  will  be  of  use  for  NDE  of  non-metallic  composite 
aerostructures. 
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2.4  Theoretical  studies 


We  have  developed  analytical  and  numerical  models  for  the  simulation  of  flaws  with  several  geometries  inside 
thick  and  thin,  current-carrying  plates  and  thin- walled  tubes.  These  models  have  been  used  for  locating  and  sizing 
flaws  in  test  samples,  and  will  serve  as  the  basis  of  future  measurement  models.  Although  the  geometry  of  real  flaws 
is  not  in  general  simple,  the  analytical  solutions  for  the  three  ideal  geometries  we  have  considered  can  be  used  for 
estimating  the  size,  orientation  and  location  of  a  flaw  faster  than  can  numerical  calculations,  and  also,  they  can  be 
used  to  determine  the  accuracy  of  more  realistic  numerical  calculations  that  describe  complicated  geometries. 


Figure  15.  SQUID  imaging  of  flaws  in  metals,  a)  calculated  for  a  9  mm  x  3  mm  ellipse  in  a  2  mm  thick 
aluminum  plate  carrying  a  current  density  of  2.5  mA/cm^.  The  major  axis  of  the  ellipse  is  horizontal,  b)  The  same 
as  (a),  except  that  the  ellipse  is  rotated  30°.  c)  Ellipse  geometry,  d)  MicroSQUID-measured  isofield  map  for  the 
ellipses  in  (c).  Adapted  from  Ref.  34. 
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2.4.1  Elliptical  holes.  We  have  calculated  two  dimensional  current  distributions  due  to  elliptical  holes  in  a 
conducting  plate,  and  have  then  used  this  distribution  to  calculate  the  magnetic  field  from  flaws  in  both  thin  and  thick 
conducting  plates.  An  analytical  solution  has  been  found  to  describe  the  magnetic  field  for  a  thick  plate'*’,  which, 
to  our  knowledge,  has  not  been  previously  reported.  This  model  can  be  used  to  determine  the  size,  location  and 
orientation  of  a  crack  in  a  current  carrying  conductor  by  setting  to  zero  the  length  of  the  minor  axis  of  the  ellipse. 


As  an  example  of  this  model.  Figs.  15a  and  b  show  the  calculated  contour  maps  of  the  magnetic  field  in  the 
z  direction  3.0  mm  above  the  plate  for  two  orientations  of  the  ellipse.  The  plate  is  2.0  mm  thick  and  is  carrying 
current  density  2.5  mA/cm^.  To  confirm  the  calculations,  we  have  mapped  the  field  for  two  ellipses  with  major  axes 
9  mm  and  12  mm,  and  minor  axes  3.0  mm,  which  are  oriented  at  90°  and  60°  with  respect  to  the  current  (see  Fig. 
15d).  While  the  ellipses  are  not  perfect  and  also  are  slightly  different  from  those  used  in  the  model  calculations,  the 
test  data  are  in  good  agreement  with  the  model  predictions.  Most  importantly,  the  model  shows  that  for  large  flaws, 
the  magnetic  field  contains  information  about  the  shape  of  the  hole  relative  to  the  direction  of  the  applied  current: 
the  fact  that  the  right-most  hole  is  tilted  with  respect  to  the  applied  current  is  evident  in  the  tilt  of  the  field  pattern. 

2.4.2  Circular  holes.  To  examine  the  effects  of  hole  size  on  the  observed  magnetic  images,  we  have  made 
an  extensive  series  of  analytical  and  finite-element  calculations  of  the  magnetic  field  produced  when  current  is  passed 
through  a  metal  sheet  containing  circular  flaws^.  We  have  previously  demonstrated  the  difference  between  imaging 
resolution  and  localizing  resolution^’,  and  we  showed  that  when  the  flaws  were  significantly  smaller  than  the  distance 
between  the  flaw  and  the  SQUID,  i.e.  the  SQUID  was  in  the  far  field  of  the  flaw,  the  magnetometer  could  localize 


Figure  16.  Proportionality  of 
magnetic  field  to  hole  area.  Varia¬ 
tion  of  the  normal  component  of  the 
magnetic  field  as  a  function  of 
hole  area  in  a  conducting  plate  at 
various  measurement  heights,  z  (in 
mm).  The  hole  radius  ranges  from 
0.1  to  3  mm.  In  the  far  field  limit 
where  the  holes  are  small  as  relative 
to  the  pickup  coil  distance  (below 
the  dashed  line),  Bj,  varies  linearly 
with  hole  area.  From  Ref.  48. 
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the  flaw  but  not  image  it.  The  most  important  finding  illustrated  in  Fig.  16  is  that  in  the  far  field  (below  the  dashed 
line),  the  peak-to-peak  amplitude  of  the  flaw  signal  was  directly  proportional  to  the  area  of  the  flaw.  In  addition, 
the  separation  between  the  positive  and  negative  extrema  in  the  magnetic  field  image  was  \f2  times  the  depth  of  the 
flaw  beneath  the  magnetometer  measurement  plane.  Thus  if  a  flaw  was  smaller  than  the  imaging  resolution  of  the 
SQUID  magnetometer,  the  flaw  size  could  still  be  determined  accurately.  For  flaws  in  the  near  field,  the  flaws  could 
actually  be  imaged  and  their  exact  shape  and  size  determined. 

2.4.3  Other  geometries.  We  have  also  developed  mathematical  models  for  simulation  of  flaws  of  different 
shapes  inside  a  thick  current-carrying  plate.  We  have  calculated  the  three-dimensional  current  distribution  due  to  a 
finite-depth  cylindrical  flaw  in  a  thick  current-carrying  plate^.  In  general,  it  is  difficult  to  solve  this  kind  of  bound¬ 
ary  value  problem  analytically.  We  have  found  an  analytic  solution  by  introducing  a  test  function,  from  which  the 
magnetic  field  is  derived,  although  the  calculations  are  rather  involved.  A  three-dimensional  current  distribution  due 
to  a  spherical  flaw  inside  a  thick  current-carrying  plate  has  been  calculated  by  the  method  of  images^'.  We  have 
found  an  analytical  solution,  which  has  not  been  reported  before,  to  calculate  the  magnetic  field  due  to  the  spherical 
flaw.  This  solution  has  been  successful  for  simulating  the  magnetic  field  produced  by  a  spherical  cavity  inside  a 
current-carrying  brass  bar,  as  shown  in  Fig.  9g. 


Tubular  samples  provide  another  challenging  geometry  which  can  be  examined  with  a  variety  of  different 
SQUID  techniques”"*'.  We  have  developed  a  lead-field  algorithm  to  calculate  the  current  distribution  from  the 
magnetic  field  produced  by  a  flaw  in  a  thin-walled,  current-carrying  cylindrical  tube*‘'“,  and  demonstrate  in  Fig.  17 
that  the  inverse  image  is  consistent  with  the  theoretical  prediction. 
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Figure  17.  Theoretical  and  experimental  images  of  current  disturbed  by  a  flaw  in  a  thin-walled  tube,  a)  The 
simulated  distribution  of  current  in  a  thin-walled  conducting  tube  containing  a  through-wall  hole  (left)  was  used  to 
determine  the  radial  component  of  the  magnetic  field  outside  the  tube  (middle),  and  a  lead  field  model  was  uaed  to 
reconstruct  the  distribution  of  current  in  the  tube  (right),  b)  The  experimental  geometry  (left),  the  measured  magnetic 
field  (middle),  and  the  image  reconstructed  from  the  data  (right).  Adapted  from  Ref.  28. 
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3.  THE  FUTURE  OF  SQUID  NDE 


3.1  The  direction  of  SQUID  NDE  research 

There  are  several  key  issues  that  need  to  be  addressed  in  SQUID  NDE.  SQUID  NDE  is  far  from  a  mature 
field:  there  is  a  vast  amount  of  work  between  the  present,  preliminary  studies,  and  a  polished,  commercial  NDE 
system  suitable  for  use  on  the  flight  line.  The  work  can  be  divided  into  at  least  three  categories:  instrumentation, 
data  acquisition  techniques  and  procedures,  and  image  processing. 

3.1.1  Instrumentation.  In  all  likelihood,  the  instrumentation  for  SQUID  NDE  will  undergo  rapid  development 
in  the  next  several  years.  Already,  compact  special-purpose  instruments  are  being  fabricated  for  use  in  industrial 
environments’’.  General  purpose  instruments  presently  under  development  will  allow  the  recording  of  susceptibility 
images  that  were  impossible  just  a  year  or  two  ago.  The  future  promises  to  provide  fully-integrated  digital  SQUID 
magnetometers,  with  self-contained  digital  feed-back  circuitry,  that  can  be  multiplexed  within  the  cryogenic  system 
and  would  require  only  minimal  room-temperature  electronics'*.  Such  SQUIDs  will  be  ideally  suited  for  imaging 
arrays,  since  an  entire  array  would  only  require  approximately  the  same  number  of  leads  as  is  presently  associated 
with  each  analog  SQUID  channel.  Smaller  SQUIDs  will  provide  even  higher  resolution  images”,  although  some 
of  this  resolution  will  be  lost  if  the  sample  cannot  be  cooled  to  the  same  temperature  as  the  SQUID.  Gradiometers 
can  be  optimized  for  specific  measurements^^””  and  can  be  configured  in  a  way  to  achieve  the  maximum 
possible  resolution  after  deconvolution  to  minimize  the  effects  of  a  finite-diameter  pickup  coil”.  Possibly  the 
greatest  short-term  improvement  will  come  not  with  isolated  improvements  to  the  SQUID  sensors  or  the  cryogenics, 
but  with  the  creation  of  optimized  systems  for  which  the  software  and  hardware  serve  as  a  seamless  and  essentially 
invisible  interface  between  the  flaw  and  the  operator  searching  for  it.  Comparison  of  today’s  laboratory  SQUID 
systems  with  commercial  eddy  current  and  ultrasound  systems  for  use  on  the  flight  line  or  on  the  oil  pipeline  should 
provide  guidance  as  to  just  how  important  the  human  factor  is  in  the  operation  of  the  instrument. 

3.1.2  Data  acquisition  techniques  and  procedures.  One  of  the  difficulties  with  SQUID  NDE  is  the  sensitivity 
of  the  system  to  fields  from  irregularities  and  discontinuities  in  the  sample  other  than  those  associated  with  the  flaw. 
The  cancellation  approach  described  above  can  successfully  reduce  the  signal  due  to  the  edge  of  the  plate;  immersing 
an  object  in  a  solution  whose  magnetic  susceptibility  matches  that  of  the  object  can  eliminate  the  susceptibility  signal 
from  the  boundary  of  an  object.  Multiple  lock-in  amplifiers  operating  at  different  frequencies  can  provide  images 
with  different  sensitivities  at  different  depths,  and  it  may  also  be  possible  to  use  an  eddy  current  approach  to 
reconstruct  the  current  distribution  as  a  function  of  depth”.  Tomographic  approaches  can  be  used  to  reconstruct 
flaws  by  applying  current  from  differing  directions.  Matched  filters  can  optimize  signals  from  expected  flaws  while 
minimizing  the  signal  from  known  inhomogeneities  characteristic  to  the  sample.  Obviously,  other  approaches  will 
be  developed  as  experience  is  gained  in  SQUID  NDE.  However,  in  assessing  the  capabilities  of  this  technique,  and 
comparing  it  to  other  approaches,  it  is  important  to  realize  that  improved  techniques  will  be  developed  not  only  for 
SQUID  NDE  but  also  other  NDE  methods.  The  hope  for  those  developing  SQUID  NDE  is  that  the  use  of  SQUIDs 
is  sufficiently  new  that  there  are  more  improvements  to  be  devised  for  SQUID  NDE  than  for  other  techniques! 

3.1.3  Image  processing.  As  we  discussed  above,  two  of  the  factors  previously  limiting  the  application  of 
SQUID  magnetometers  to  the  nondestructive  testing  of  materials  and  structures  have  been  the  lack  of  sufficient  spatial 
resolution  and  techniques  to  provide  images  of  flaws.  In  the  past  several  years,  SQUID  NDE  has  made  great  progress 
in  overcoming  these  limitations.  However,  due  to  the  physical  configuration  of  SQUID  magnetometers  that  measure 
the  field  from  warm  samples,  there  is  an  inherent  limitation  in  SQUID  spatial  resolution  capability.  In  addition,  as 
with  any  real  world  measurement  procedure,  the  SQUID  measurement  process  introduces  other  forms  of  imaging 
distortion.  In  order  to  remove  this  distortion  and  to  increase  resolution,  it  is  possible  to  apply  the  concept  of  signal 
restoration  to  the  measured  images.  The  objective  of  a  signal  restoration  algorithm  is  to  reconstruct  (or  recover)  an 
image  from  a  degraded  measure  of  that  image  by  employing  a  priori  knowledge  of  the  distortion  process.  A  critical 
step  in  any  useful  signal  restoration  algorithm  is  that  of  generating  a  reliable  model  of  the  measurement  distortion 
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process.  A  restored  image  can  be  obtained  in  principal  by  applying  an  inverse  of  this  distortion  model  to  the 
measured  image.  Our  comparisons  between  numerical  models  and  experimental  data  suggest  that  SQUID  NDE  has 
the  great  advantage  of  providing  reproducible,  quantitative  data  well  suited  for  image  restoration  and  quantification 
of  flaw  size.  The  challenge  is  to  optimize  the  analysis  and  flaw  identification  process. 

Many  of  the  earlier  forms  of  image  restoration  were  carried  out  in  the  spatial  frequency  domain^’,  as  was 
done  for  Fig.  2.  The  two-dimensional  Fourier  transform  normally  played  an  important  role  in  such  restoration  algo¬ 
rithms.  However,  Fourier-transform-based  techniques  have  an  inherent  resolution  limitation  which  may  limit  their 
use  in  applications  requiring  high  resolution.  Furthermore,  it  is  difficult  to  apply  a  priori  constraints  to  the  image, 
such  as  the  requirement  that  no  current  extend  beyond  the  known  borders  of  the  sample.  It  is  for  this  reason  that 
the  spatial  domain  became  the  natural  terrain  for  many  of  the  more  recently  developed  image  restoration  algorithms. 
In  a  spatial  domain  approach,  as  in  the  finite  element  reconstruction  shown  in  Fig.  3,  the  need  to  properly  model  the 
distortion  process  remains  an  essential  ingredient.  A  typical  spatial-domain-based  image  restoration  algorithm  is 
heavily  based  on  modern  linear  algebraic  techniques,  and  possesses  the  ability  to  effect  "super  resolution"  image 
restorations’®.  Undoubtedly,  there  will  be  trade-offs  dictated  by  measurement  noise,  but  it  is  reasonable  to  assume 
that  such  advanced  algorithms  will  be  able  to  provide  noticeable  enhancement  of  the  images  over  what  we  have  shown 
so  far. 

3.2  SQUID  NDE  versus  other  techniques 

The  second  key  concern  is  to  further  our  understanding  of  the  capabilities  and  limitations  of  the  technique 
relative  to  other  NDE  techniques.  This  will  require  detailed,  quantitative  comparisons  of  results  obtained  on  a 
common  set  of  samples.  Such  an  analysis  of  SQUID  NDE  could  be  aided  by  a  theoretical  analysis  of  the  sensitivity 
of  a  SQUID  to  a  particular  type  of  flaw  measured  in  a  certain  manner.  The  development  of  numerical  measurement 
models  that  include  not  only  SQUID  noise,  but  also  environmental  noise  and  realistic  sample  inhomogeneities  and 
geometries  will  prove  to  be  challenging,  but  will  provide  a  valuable  estimate  of  the  smallest  flaw  that  could  be 
detected  and  the  largest  flaw  that  might  be  missed.  The  actual  determination  of  a  probability  of  detection  (POD)  for 
SQUID  NDE  of  a  given  class  of  flaws  will  require  a  substantial  experimental  effort.  Although  one  can  always  hope 
to  identify  a  particular  flaw  that  can  be  detected  with  SQUIDs  and  not  by  any  other  competing  technique,  the  general 
acceptance  of  SQUIDs  as  regular  tools  for  aircraft  NDE  will  more  likely  be  contingent  on  whether  such  a  careful 
POD  analysis  proves  SQUID  NDE  to  be  better  than  competing  techniques,  at  least  for  a  subset  of  known  aircraft 
flaws. 


If  it  can  be  shown  in  the  laboratory  that  SQUIDs  have  a  potential  role  to  play  in  aircraft  NDE,  we  must  then 
address  the  constraints  imposed  by  the  cryogenic  requirements  and  cost  of  SQUID  systems.  At  present,  SQUIDs  are 
harder  to  use  than  eddy  current  techniques:  the  probes  are  larger,  heavier,  and  require  scanning  systems.  They  are 
exquisitely  sensitive  to  magnetic  noise  from  DC  to  radio  frequencies.  The  systems  are  presently  expensive,  costing 
between  $25,000  and  $500,000  each,  depending  upon  the  complexity  and  capabilities  of  the  system.  SQUID  NDE 
systems  are  presently  comparable  in  cost  and  size  to  ultrasound  C-scan  systems  with  water  baths  and  X-ray 
tomography  units.  The  use  of  liquid  helium,  while  of  little  concern  to  a  physicist  in  a  research  laboratory,  might 
prove  to  be  an  impediment  to  use  of  the  instrument  at  an  Air  Force  Air  Logistic  Center  or  a  commercial  airline 
maintenance  depot.  The  interpretation  of  the  data  at  present  requires  sophisticated  mathematical  models.  These  are 
not,  however,  insurmountable  problems,  but  instead  challenges  to  the  clever  systems  engineer.  With  the  advent  of 
high  transition-temperature  superconductivity  (HTS),  SQUIDs  and  their  simplified  cryogenic  apparatus  will  quickly 
become  smaller,  lighter,  and  less  expensive.  Because  of  the  large  amplitude  of  many  NDE  signals,  the  somewhat 
lower  sensitivity  of  HTS  SQUIDs  as  compared  to  helium-temperature  ones  is  not  a  major  issue.  Multiple  channel 
SQUIDs,  gradiometers  designed  to  detect  flaws  but  not  background  signals,  and  user-friendly  electronics  might 
eventually  make  SQUIDs  as  easy  to  use  as  the  Magneto-Optical  Inspection  technique,  but  vastly  more  sensitive.  It 
is  also  important  to  recognize  that,  in  contrast  with  established  NDE  technologies,  the  present  cost  of  SQUID  systems 
does  not  reflect  any  economies  of  scale,  since  most  systems  sold  so  far  have  been  custom-built  units. 
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While  general-purpose  SQUID  systems  might  be  useful,  the  tradeoffs  of  cost,  reliability,  and  ease  of  use 
suggest  that  the  special-purpose  SQUID  optimized  for  a  single  type  of  flaw  and  a  particular  measurement  may  lead 
to  more  rapid  introduction  of  SQUIDs  onto  the  flight  line.  In  that  case,  we  need  to  identify  the  type  of  flaw  for 
which  the  early  flight-line  prototypes  should  be  developed. 

Unfortunately,  once  the  technical,  cost,  and  practical  limitations  of  SQUIDs  are  overcome,  the  acceptance 
of  this  technology  by  the  aircraft  user  and  maintenance  communities  is  not  guaranteed.  The  aircraft  community  is 
historically  conservative  with  regard  to  NDE  technologies,  driven  in  part  by  concerns  that  the  premature  adoption 
of  a  new  technology  might  lead  to  litigation  should  an  undetected  flaw  contribute  to  an  accident.  If,  however,  a 
SQUID  NDE  system  can  be  developed  that  is  simple  to  use,  inexpensive,  is  possibly  optimized  for  a  single  type  of 
flaw,  and  can  detect  that  type  of  flaw  with  unprecedented  sensitivity  and  reliability,  then  the  acceptance  of  SQUID 
NDE  will  be  much  more  likely. 

3.3  Promising  applications 

Given  what  we  have  already  learned  and  the  expected  improvements  in  SQUID  NDE  instrumentation  and 
techniques,  it  is  reasonable  to  hypothesize  as  to  the  areas  in  which  SQUID  NDE  might  have  the  greatest  impact.  At 
present  we  see  several,  including  second-layer  cracks  and  corrosion  in  aluminum  fuselage  and  wing  coverings  and 
susceptibility  imaging  of  composites  and  magnetic  decorants. 

3.3.1  Second  layer  cracks  and  corrosion.  The  preliminary  data  presented  above  supports  our  hypothesis  that 
low-frequency  eddy-current  NDE  with  SQUIDs  is  well  suited  for  detecting  second-layer  cracks.  Other  techniques 
have  sufficient  limitations  to  justify  much  more  research  in  this  area:  ultrasound  is  capable  of  detecting  flaws  in 
second  and  third  layers  only  if  there  is  an  intact  layer  of  bonding  material  between  each  layer;  thermography  is  well 
suited  for  detecting  disbonds  between  these  layers,  but  not  cracks  at  rivets;  the  MOI  system  cannot  at  present  operate 
at  low  enough  frequencies  to  see  deep  into  metallic  structures. 

3.3.2  Measurements  of  active  corrosion.  SQUID  magnetometry  has  seen  little  application  to  the  study  of 
active  corrosion  in  metals,  primarily  because  of  the  low  spatial  resolution  provided  by  SQUIDs  optimized  for  bio- 
magnetic  measurements  and  the  unfamiliarity  of  corrosion  scientists  with  SQUIDs.  The  most  definitive  SQUID  study 
to  date  of  corrosion  was  performed  by  Bellingham^’,  and  demonstrated  the  ability  of  a  low-resolution  SQUID  to  map 
corrosion  currents  in  a  dish.  More  substantive  measurements  utilized  the  SQUID  magnetometer  solely  as  a  low- 
impedance  ammeter  and  examined  the  factors  that  affected  amplitude  and  spectrum  of  corrosion  noise  in  a  zinc-hydro¬ 
chloric  acid  nonvoltaic  cell.  These  particular  measurements  might  have  been  made  more  easily  using  a  low  imped¬ 
ance,  low  noise  current-to-voltage  converter*”.  High  resolution  studies  of  corrosion  might  be  facilitated  by  the  use 
of  a  high-resolution  SQUID  magnetometer,  such  as  MicroSQUID^  or  the  system  recently  manufactured  by  Quantum 
Magnetics*' . 

The  ability  of  SQUID  magnetometers  to  detect  steady  magnetic  fields  makes  them  ideal  for  mapping  of  inter¬ 
nal  corrosion  currents.  Mathematical  models  developed  by  the  Vanderbilt  group  to  describe  current  flow  in  aniso¬ 
tropic  cardiac  tissue*^,*^  are  in  fact  applicable  to  the  distribution  of  currents  in  advanced  metallic  composite  conduc¬ 
tors.  With  a  high  resolution  SQUID  such  as  MicroSQUID,  it  would  be  possible  to  scan  an  otherwise  non-magnetic 
composite  sample  of  either  regular  or  irregular  shape,  to  determine  whether  steady  magnetic  fields  are  produced  by 
corrosion.  A  casual  experiment  at  Vanderbilt  demonstrated  that  MicroSQUID  was  able  to  detect  large  DC  fields  pro¬ 
duced  by  two  moistened  layers  of  aluminum  spanned  by  two  dissimilar  metal  rivets.  The  studies  by  Bellingham  and 
MacVicar*^  indicate  that  SQUID  magnetometers  also  have  the  adequate  sensitivity  and  spatial  resolution  to  measure 
the  frequency  spectrum  of  corrosion  currents  that  arise  between  a  metal  and  a  strong  electrolyte;  a  preliminary  study 
at  Vanderbilt  demonstrates  that  MicroSQUID  can  image  the  even  weaker  Johnson  noise  arising  from  thermal  motion 
of  electrons  in  a  ring  of  copper.  While  it  is  unlikely  that  Johnson  noise  images,  such  as  that  shown  in  Fig.  18,  will 
be  of  practical  utility  for  flight  line  NDE  inspections,  these  data  clearly  demonstrate  that  SQUIDs  have  adequate  sensi- 
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tivity,  spatial  resolution,  and  frequency  response  to  image  both  the  DC  component  and  the  full  AC  spectrum  of 
currents  from  ongoing  corrosion. 


Because  the  magnetic  fields  from  such  corrosion  currents  are  weak  compared  with  the  earth’s  magnetic  field, 
and  may  be  subject  to  contamination  by  thermoelectric  currents  in  entire  aircraft  structures,  it  is  likely  that  SQUID 
techniques  might  never  be  useful  for  measuring  corrosion  currents  on  the  flight  line.  However,  this  technique  pro¬ 
vides  a  new  and  potentially  powerful  tool  for  studying  corrosion  of  aircraft  samples  in  the  laboratory.  The  ability 
to  measure  the  strength  and  spatial  distribution  of  filiform  currents  beneath  a  layer  of  paint,  and  corrosion  currents 
from  nonvisible  corrosion  between  the  first  and  second  or  second  and  third  layers  of  a  sample  would  provide  a  new, 
noninvasive  monitoring  technique  for  evaluating  the  effectiveness  of  corrosion  prevention  strategies  and  for  mon¬ 
itoring  corrosion  development  in  both  time  and  space. 


Figure  18.  Images  of  Johnson  Noise.  Images  of  the  magnetic  field  produced  by  Johnson  noise  currents  flowing 
in  a  copper  ring,  a)  Images  for  the  0  to  6  Hz  band,  b)  Images  for  the  4994  to  5000  Hz  band.  Note  from  the 
magnetic  field  distribution  that  the  low  frequency  currents  flow  around  the  ring  but  the  high  frequency  ones  only  flow 
locally. 


3.3.3  Susceptibility  Imagine.  The  ability  of  magnetic  susceptibility  tomography  to  distinguish  materials  on 
the  basis  of  the  diamagnetic  and  paramagnetic  susceptibilities  may  allow  for  discrimination  of  materials  based  upon 
their  chemical  makeup.  For  example,  the  magnetic  susceptibility  of  aluminum  is  -f-20.9  x  10^  (Sl/m^),  while  that 
of  aluminum  oxide  is  -18.1  x  10*,  producing  a  susceptibility  contrast  four  times  greater  than  the  readily  imaged 
susceptibility  of  plexiglas.  For  planar  samples,  scanning  susceptometry  would  provide  a  real-time  2-D  susceptibility 
image.  For  three-dimensional  objects,  susceptibility  tomography  should  be  able  to  identify  corrosion  products  deep 
inside  the  sample.  The  ability  to  image  small  quantities  of  ferromagnetic  or  superparamagnetic  tracers  offers  promise 
for  rapid  detection  of  surface-breaking  flaws  in  ceramic  and  other  non-metallic  materials,  as  well  as  in  non¬ 
ferromagnetic  metals.  Since  MicroSQUID  with  room-temperature  magnetizing  coils  does  not  represent  an  optimized 
imaging  susceptometer,  the  development  of  this  technique  will  require  a  totally-superconducting  susceptometer.  The 
increased  sensitivity  of  such  a  system  should  provide  previously  unimagined  capabilities  for  SQUID  NDE. 
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The  future  of  SQUID  NDE  for  aircraft  is  still  uncertain,  largely  because  there  have  yet  to  be  any  proven 
applications  for  which  SQUID  NDE  is  clearly  better  than  competing  techniques.  But  judging  from  the  progress  that 
has  been  made  in  the  past  several  years  and  the  potential  for  further  improvements,  SQUID  NDE  may  make  its  appear¬ 
ance  on  the  flight  line  within  the  next  few  years.  The  extent  to  which  SQUID  NDE  can  make  real  inroads  into  the 
conventional  NDE  market  will  depend  in  part  on  the  success  of  the  ongoing,  initial  attempts  at  transferring  this 
technology  out  of  the  laboratory.  The  research  reviewed  in  this  paper  should  provide  some  insights  as  to  the  direction 
the  field  has  been  taking,  and  a  rather  speculative  view  of  the  future.  All  in  all,  the  future  of  SQUID  NDE  for  aircraft 
seems  bright,  but  there  is  much  work  to  be  done  before  this  can  be  assured. 
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